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Abstract 

The computation of long wave propagation through the ocean obvi- 
ously depends on the initial condition. When the waves are generated 
by a moving bottom, a traditional approach consists in translating the 
"frozen" sea bed deformation to the free surface and propagating it. 
The present study shows the differences between the classical approach 
(passive generation) and the active generation where the bottom mo- 
tion is included. The analytical solutions presented here exhibit some 
of the drawbacks of passive generation. The linearized solutions seem 
to be sufficient to consider the generation of water waves by a moving 
bottom. 

Resume 

Theorie lineaire de generation de vagues par mouvement du 

fond. Les calculs de propagation d'ondes longues a travers I'ocean 
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doivent naturellement etre alimentes par la condition initiale. Le but 
de cette note est de montrer I'insufRsance de I'approche classique qui 
consiste a translate! la deformation "gelee" du fond vers la surface libre 
et a la laisser se propager. Un calcul analytique qui met en evidence 
les inconvenients de I'approche classique de generation passive est 
presente ici. Les solutions linearisees semblent etre bien adaptees pour 
traiter la generation de vagues par mouvement du fond. 

Mots-cles : Vagues linearisees, Mecanique des fluides, Probleme de 
Cauchy-Poisson, Generation des tsunamis 



Version francaise abregee 

Le probleme de la generation des tsunamis est un sujet relativement 
recent. L'un des pionniers dans ce domaine fut Hammack [1]. Le but de 
cette note est d'apporter une contribution a ce probleme. La condition ini- 
tiale utilisee dans les codes de propagation des tsunamis est souvent obtenue 
en translatant a la surface libre la deformation du fond suite a un tremblement 
de terre. Cette approche presente plusieurs inconvenients. Tout d'abord la 
dynamique du processus de generation est negligee. II est evident qu'un 
glissement lent ne produit pas des vagues de meme amplitude qu'un glisse- 
ment rapide. Ensuite le champ des vitesses initiales est egalement neglige. 
Nous presentons dans cette note un simple modele de generation de tsunamis 
(voir [3] pour plus de details). 

Le probleme des ondes de surface est tout d'abord linearise (l)-(3) pour 
un fond qui a un mouvement prescrit (generation active). II est resolu par 
la methode des transformees de Laplace en temps et de Fourier en espace. 
En supposant que le mouvement du fond est instantane, on obtient pour la 
deformation de la surface libre rii{x,y,t) I'expression (10). Les vitesses peu- 
vent egalement etre calculees. Le probleme des ondes de surface est ensuite 
linearise dans le cas oil la deformation du fond est simplement translatee 
jusqu'a la surface libre (generation passive). Ce probleme est different du 
precedent. En effet, puisque le fond reste immobile en tout temps, la condi- 
tion cinematique au fond devient (11). Par ailleurs, la condition initiale sur la 
surface libre devient r]{x, y, 0) = ({x, y), oii ({x, y) represente la deformation 
permanente du fond. On obtient alors pour la deformation de la surface libre 
r]{x,y,t) I'expression (16). 

On compare ensuite les deux solutions (10) et (16). Dans les deux cas on 
utilise la meme deformation du fond due a un tremblement de terre, qui est 
donnee par la solution d'Okada [2] pour une faille finie rectangulaire situee 
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a une profondeur de 3 km, de longueur 6 km et de largeur 4 km. Les autres 
parametres sont: module d'Young = 9.5 GPa, coefficient de Poisson = 0.27, 
glissement = 15 m. La profondeur d'eau est 1 km et I'acceleration de la 
gravite est g = 10 m/s^. Pour ce type de deformation du fond, la vague 
initiale a la forme d'un et I'axe y represente la direction privilegiee pour 
la propagation des ondes. La figure (1) montre le proffi des vagues mesure 
a plusieurs endroits le long de la surface libre. La courbe en trait plein 
represente la solution dynamique tandis que la courbe en traits pointilles 
represente le scenario de generation passive. Les amplitudes sont clairement 
plus grandes dans ce dernier C3,S. Let figure (2) montre la difference relative 
(17) entre les deux solutions. II y a deux differences essentielles entre les 
deux solutions: la generation passive donne des amplitudes de vagues plus 
elevees et dans le cas de la generation active la colonne d'eau joue le role 
d'un filtre qui attenue les hautes frequences grace a la presence du cosinus 
hyperbolique au denominateur. Les resultats dependent naturellement de 
I'echelle de temps caracteristique de la deformation du fond. Dans le futur, 
nous allons egalement etudier I'effet des termes non-lineaires. 

1 Introduction 

Tsunami generation is a relatively recent topic inspired for example by 
the pioneer work of Hammack [1]. Since then, progress has been moderate. 
The present note provides a contribution to the development of this field of 
hydro dynamics . 

The computation of long wave propagation across the ocean is a com- 
plicated task. The accuracy of the results depends on different factors such 
as the numerical method, the discretization error, the mathematical model 
error and others. The error made in the initial condition cannot be corrected 
by the numerical method and will propagate in space and time. In our opin- 
ion it is important to construct an initial condition that is as accurate as 
possible. Surprisingly there has been relatively little research in this field. 

The initial condition is often constructed as follows. One takes coseismic 
deformations predicted by various models ([2] is presently used in many cases) 
and translates them to the free surface. The velocity field is assumed to be 
zero. Then, a finite difference code computes the gravitational wave train 
induced by this free-surface disturbance. 

This approach has several drawbacks. First of all, the dynamic character 
of the tsunami generation process is not taken into account. It is obvious 
from physical intuition (and con&med by relatively simple computations) 
that slow slip does not produce waves of the same amplitude as fast bottom 
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motion. So, usually, the initial wave amplitude is either under- or over- 
estimated depending on the time characteristics of the source. Moreover 
the initial velocity field in the fiuid due to the moving bottom is completely 
neglected. Our computations show that this is not necessarily true. 

The present note sheds some light on these drawbacks. At the same time 
the model studied in the present note can be considered as one of the simplest 
models for dynamic tsunami generation. We refer to [3] for more details. 

2 Linearized waves 

2.1 Moving bottom solution 

Let us consider a three-dimensional fiuid domain Q bounded above by 
the free surface of the ocean z = ri{x, y, t) and below by the rigid ocean fioor 
z = —h + ({x, y, t). The domain Q is unbounded in the horizontal directions 
X and y, and can be written as f2 = x [—h + C{x,y,t),ri{x,y,t)]. It is 
assumed that the fiuid is incompressible and the fiow irrotational. The latter 
implies the existence of a velocity potential y, z, t) which completely 
describes this fiow. Initially the fiuid is assumed to be at rest and the sea 
bottom to be horizontal [z = —h). Mathematically these conditions can 
be written in the form of initial conditions ^, 0) = 0, ri{x,y,0) = 

and ({x,y,0) = 0,^ which complete the formulation of the initial boundary 
value Cauchy-Poisson problem described below. Thus, at time t = 0, the free 
surface and the sea bottom are defined by 2; = and z = —h, respectively. 
At time t > the bottom boundary moves in a prescribed manner which is 
given by 2; = —h + ({x, y, t). The displacement of the sea bottom is assumed 
to have all the properties required to compute its Fourier transform in x, y 
and its Laplace transform in t. The resulting deformation of the free surface 
z = ri{x, y, t) must be found. 

Solving this problem is a difficult task due to the nonlinearities and the a 
priori unknown free surface. In this study we linearize the equations and the 
boundary conditions. The linearized problem in dimensional variables reads 



A0 = 



ix,y,z)eR'x[ 



h,0], 



(1) 



(90 drj 

dz dt ' 




z = 0, 



(2) 



^Thc last condition is not an initial condition. Wc added it in order to have a flat 
bottom initially. In fact, it is not required for the mathematical method. 
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Combining equations (2) yields the single free-surface condition 

d'^d) d(h , , 

g| + 9a^ = o. . = 0. (4) 

The problem (1), (3), (4) can be solved by using the method of inte- 
gral transforms. We apply the Fourier transform in (x, y) and the Laplace 
transform in time t. For the combined Fourier and Laplace transforms, the 
notation F{k, i, s) is introduced. After applying the transforms, equations 
(1), (3) and (4) become 

^-(P + a^ = 0, (5) 



^(k,i,-h,s) = sak,i,s), (6) 

s^^{kj,0,s) + g^{kj,0,s) = 0. (7) 
The transformed free-surface elevation can be obtained from (2): 

rj{k,i,s) = --^ik,i,0,s). (8) 
9 

A general solution of equation (5) is given by 

£, z, s) = A{k, £, s) cosh(m2;) + B{k, i, s) sinh(m2;), (9) 



where m = The functions A{k,i,s) and B{k,i,s) can be easily 

found from the boundary conditions (6) and (7): 



cosh(m/i) [s^ + gm tanh(m/i)] ' 



B{k,i,s) 



m cosh(m/i) [s^ + gm tanh(m/i)] 



From now on, the notation u = ^^/ gmtanh.{mh) will be used. Substituting 
the expressions for the functions A, B in the general solution (9) yields 

(p{k,e,z,s) = — - cosh{mz) smh(m2;) . 

cosh(m/i)(s"' + u;^) \ gm J 
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From (8), the free-surface elevation becomes f]{k,£,s) = s'^C{k, i, s){s'^ + 
uj'^)^^/ cosh(m/i). 

Now we assume that the sea bed deformation is instantaneous, i.e. ({x, y, t) 
C{x, y)H(t), where H(t) denotes the Heaviside step function^. After some an- 
alytic computations one obtains the final integral formula for the free-surface 
elevation: 

r]i{x,y,t) = —— // \ cos ujt dkdi, (10) 

(zvr)"^ J J cosh(m/i) 

where ({k, i) is the Fourier transform of ({x, y). The velocity field due to the 
moving bottom can also be computed [3]. It cannot necessarily be neglected. 

2.2 Passive generation 

In this case the initial condition is obtained by translating the sea bed 
deformation to the free-surface and the evolution of this system is computed. 
Next we give an analytic solution to this problem. This solution is supposed 
to model what happens in the classical tsunami generation approach. 

First of all, we have to make several modifications to the previous prob- 
lem. Since the sea bed remains fixed at all time, the kinematic condition at 
the bottom becomes 

l^^o, ,U) 

The main difference with §2.1 is the initial condition on free surface which 
becomes ri{x, y, 0) = ({x, y). 

Again we apply the Fourier transform in the horizontal coordinates with 
the notation F{k, €). We do not apply the Laplace transform because there is 
no substantial dynamics in this problem. Equation (5) is the same as before 
while (7) and (11) become 



a. 



g^^ik,i,0,t) + g-^{k,i,0,t)=0, (12) 

1^ = 0, z = -h. (13) 

Since Laplace's equation still holds we have the same general solution 
(9). The relation between the functions A{k,i,t) and B{k,i,t) can be easily 



^Tlie Heaviside funetion has the property that it is equal to zero for t < 0. So, choosing 
this particular form for C satisfies automatically the condition ({x,y,0) = 0. 
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found from the boundary condition (13): 

B{k, i, t) = A{k, i, t) tanli(m/i). (14) 
From equation (12) and the initial conditions one finds A{k,i,t) so that 

(f){k,i,z,t) = ——({k,£) smujt(^cosh.{mz) + tanh(m/i) sinh(m2;) j . (15) 

From the transformed dynamic condition 0j + (yfr/ = at 2; = 0, it is easy 
to find the Fourier transform of the free surface elevation 

rj{k, i, t) = C{k, i) cos out. 

The inversion of the Fourier transform provides the simple integral solu- 
tion 

V{x,y,t) = j j ak,i)cosujte'^''''+'y^dkd£. (16) 

2.3 Numerical computation 

We now compare the two solutions (10) and (16). In both cases we use the 
same sea bed deformation due to an earthquake which is given by Okada's 
solution [2] for a finite rectangular fault occurring at depth 3 km with length 
6 km and width 4 km. The other parameters are: Young's modulus = 9.5 
GPa, Poisson's ratio = 0.27, dip angle = 13°, strike angle = 90° and slip 
= 15 m (dip-slip faulting). The water depth is 1 km and the acceleration 
due to gravity is = 10 m/s^. For this particular sea bed deformation, the 
initial wave has a A^— shape and the |/— axis is the preferred direction for 
wave propagation. 

All integrals were computed with a Filon-type numerical integration for- 
mula [4], which takes into account the oscillatory behaviour of the integrands. 

Figure (1) shows the wave profile measured at several locations along the 
free surface. The solid line represents the solution with instantaneous bottom 
deformation while the dashed line represents the passive wave generation 
scenario. The latter clearly exhibits higher wave amplitudes. 

Figure (2) represents the relative difference between the two solutions 
which is defined by 

/ .^ \Viix,y,t) -rj{x,y,t)\ 
r{x,y,t) = r—r- . (17) 

I WiWoo 

Intuitively this quantity represents the deviation of the passive solution from 
that generated by a moving bottom in units of the maximum amplitude of 
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Figure 1: Wave profile ri{x,y,t) in km along several artificial tide gauges 
(x, y) in km versus time t in s 
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Figure 2: Relative difference between the solutions (10) and (16) correspond- 
ing to active and passive generation versus time 
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3 Conclusions 

Looking at the analytic expressions for r] and the numerical results gives 
some clear conclusions. Let us focus on two main differences which can be 
crucial for accurate tsunami modelling. 

First of all, the wave amplitudes obtained with the instantly moving bot- 
tom are lower than those generated by initial translation of the bottom mo- 
tion (this statement follows from the inequality cosh mh > 1 and a compari- 
son between formulas (16) and (10)). The numerical experiment shows that 
this difference is typically of the order of 20%. 

The second feature is more subtle. The water column has an effect of low- 
pass filter. It means that if the initial deformation contains high frequencies 
they will be attenuated in the moving bottom solution because of the hyper- 
bolic cosine cosh(/ivP+^) in the denominator that grows exponentially 
with m. 

Let us mention that if we prescribe a more realistic bottom motion as 
in [3] for instance, the results will depend on the characteristic time of the 
sea-bed deformation. Even for very fast bottom motions, the generated wave 
amplitude will never reach the passive generation solution. For slow motions, 
the amplitude will be in general much smaller. 

Future studies will provide a more thorough development of this topic 
including the effect of different nonlinearities. 
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